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The momentum distributions, natural orbits, spectroscopic factors and quasi-hole wave functions 
of the ^^C ,^®0 ,*°Ca ,^*Ca , and ^"^Pb doubly closed shell nuclei, have been calculated in the 
framework of the Correlated Basis Function theory, by using the Fermi hypernetted chain resumma- 
tion techniques. The calculations have been done by using the realistic Argonne Wg nucleon-nucleon 
potential, together with the Urbana IX three-body interaction. Operator dependent correlations, 
I which consider channels up to the tensor ones, have been used. We found noticeable effects pro- 

. duced by the correlations. For high momentum values, the momentum distributions show large 

' enhancements with respect to the independent particle model results. Natural orbits occupation 

04 ' numbers are depleted by about the 10% with respect to the independent particle model values. The 

5^ , effects of the correlations on the spectroscopic factors are larger on the more deeply bound states. 
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y—i ^ One of the major achievements of nuclear structure studies in the last ten years is the consolidation of the validity 
of the non relativistic many-body approach. The idea is to describe the nucleus with a Hamiltonian of the type: 

o ■ 

o: ^ = -^E^'+ E ^«^+ E w 

' i i<j = l i<j<k=l 

o ■ 

' where the two- and three-body interactions, Vij and Vijk respectively, are fixed to reproduce the properties of the 
. two- and three-body nuclear systems. The Schrodinger equation has been solved without approximations for few 
,_L( ' body systems [l| and light nuclei 0, up to A=12 The obtained results provide good descriptions, not only of the 
O [ energies of these nuclei, but also of other observables. 

The difficulties in extending to medium and heavy nuclei the techniques used in few body systems and light nuclei, 
favored the development of models, and of effective theories. The basic idea of the effective theories is to work in 
a restricted space of the many-body wave functions. Usually, one works with many-body wave functions which are 
Slater determinants of single particle wave functions. The idea of single particle wave functions implies the hypothesis 
of a mean-field where each nucleon move independently from the other ones. This Independent Particle Model (IPM) 
is quite far from the picture outlined by the microscopic calculations quoted above, describing the nucleus as a many- 
body system of strongly interacting nucleons. In the Hartree-Fock theory, which provides the microscopic foundation 
of the IPM, the Hamiltonian is not any more that of Eq. ([1]), but it is an effective Hamiltonian built to take into 
account, obviously in an effective manner, the many-body effects that the microscopic calculations explicitly consider. 
The construction of effective interactions starting from the microscopic ones, covers a wide page of the nuclear physics 
history, starting from the Brueckner G-matrix effective interactions [1, @, up to the recent Viowk interaction [1,01 and 
the interaction obtained by applying the unitary correlation operator method [1, . 

The application of the IPM is quite successful, but there are evidences of the intrinsic limitations in its applicability. 
For example, the measured spectroscopic factors are systematically smaller than one [13, EH 113 > which is the value 
predicted by the IPM. The (e,e'p) cross sections in the quasi-elastic region need a consistent reduction of the IPM 
hole strength to be explained [ij, [l3| ■ The same holds for the electromagnetic form factors of the low- lying states, 
especially those having large angular momentum 1^ 



16 L The emission of two like nucleons in photon and electron 



scattering process cannot be described by the IPM [IT], [18| . Also the charge density distributions extracted by elastic 
electron scattering data are, in the nuclear interior, smaller than those predicted by the IPM [19, 20]. These examples 
indicate the presence of physics effects, commonly called correlations, which are not described by the IPM. 

It is common practice to distinguish between long- and short-range correlations since they have different physical 
sources. The long-range correlations are related to collective excitations of the system, such as the giant resonances. 
The short-range correlations (SRC) are instead connected to the strongly repulsive core of the microscopic nucleon- 
nucleon interaction. The repulsive core reduces the possibility that two nucleons can approach each other, and this 
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modifies the IPM picture where, by definition, the motion of each nucleon does not depend on the presence of the 
other ones. 

Even though most of the calculations of medium heavy nuclei are based on the IPM, and on the effective theories, 
various techniques, aiming to attack the problem by using the microscopic Hamiltonian ([1]), have been developed. The 
Brueckner-Hartree-Fock approach has been recently applied to the ^^O nucleus [21]. No core-shell model calculations 
have been done for nuclei lighter than ^^C [12, [2^. The coupled cluster method has been used to evaluate ^^O 
properties [13, [111 . 

About fifteen years ago |26j, we started a project aimed to apply to the description of medium and heavy nuclei 
the Correlated Basis Function (CBF) theory, successfully used to describe the nuclear and neutron matter properties 
[13, HI]. We solve the many-body Schrodinger equation by using the variational principle: 

SE[^] = ' , ' ^ = . (2) 

The search for the minimum of the energy functional is done within a subspace of the full Hilbert space spanned by 
the A-body wave functions which can be expressed as: 

^{A)^J^{1,...,AM1,...,A), (3) 

where A) is a many-body correlation operator and is a Slater determinant composed by single 

particle wave functions, (j)a{i)- In our calculations, we used two-body interactions of Argonne and Urbana type, 
and we considered all the interaction channels up to the spin-orbit ones. Together with these two-body interactions, 
we used the appropriated three-body forces of Urbana type. The complexity of the interaction required the use of 
operator dependent correlations. We consider correlations of the type: 

^-s[ n F,A , (4) 

where 5 is a symmetrizer operator and Fij is expressed in terms of two-body correlation functions fp as: 

6 

In the above equation we have adopted the nomenclature commonly used in this field, by defining the operators as: 

0^='''^[l,cT, -a,,S,,]^[l,T,-T,] , (6) 

where cr^ and indicate the usual Pauli spin and isospin operators, and Sij — (3fij • CTify • (Tj — cr^ • crj) is the tensor 
operator. 

We recently succeeded in formulating the Fermi Hypernetted Chain (FHNC) equations, in Single Operator Chain 
(SOC) approximation, for nuclei non saturated in isospin, and with single particle basis described in a jj coupling 
scheme. We presented in Ref. ^29.] the binding energies and the charge distributions of ^^C , ^^O , '"'Ca , "^^Ca and 
208pj^ doubly closed shell nuclei obtained by using the minimization procedure These calculations have the same 
accuracy of the best variational calculations done in nuclear and neutron matter [2?! . [28j . 

In the present article, we show the results of our study, done in the FHNC/SOC computational scheme, on some 
ground state quantities related to observables. They are momentum distributions, natural orbits and their occupation 
numbers, quasi-hole wave functions and spectroscopic factors. We used the many-body wave functions obtained in 
Ref. by solving Eq. ^ with the Ar gonne Wg two-nucleon potential, together with the Urbana IX three-body 
force. We have calculated momentum distributions also with the wave functions produced by another interaction, 
the Urbana U14 truncated up to the spin-orbit terms, implemented with the Urbana VII three-body force. The 
results obtained with this last interaction do not show relevant differences with those obtained with the Wg and UIX 
interaction, therefore we do not present them. 

The paper is organized as follows. In Sect. [n]we present the results of the One-Body Density Matrix (OBDM) 
and of the momentum distribution. In Sect. IIIII we discuss the natural orbits, i.e. the single particle wave functions 
forming the basis where the OBDM is diagonal. In Sect. II VI we present our results about the quasi-hole wave functions 
and in Sect. |V]we summarize our results and draw our conclusions. 
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II. ONE-BODY DENSITY MATRIX AND MOMENTUM DISTRIBUTION 

We define the one-body density matrix, (OBDM), of a system of A nucleons as: 
pixux[) ^ ^p^'^'^*(ri,rl)xt(l)xl(l)x.'(l')x.(l') 

s.s' .t 

= ^ ,^|,^ ^ j dx2...dXA'^\xi,X2,...,XA)^(x'-i^,X2,---,XA) ■ (J) 

In the above expression, the variable Xi indicates the position (r^) and the third components of the spin (s^) and 
of the isospin (ti) of the single nucleon. The x(i) functions represent the Pauli spinors. With the integral sign we 
understand that also the sum on spin and isospin third components of all the particles from 2 up to A, is done. In 
our calculations we are interested in the quantity: 

p\r,,r[)= [p^^-\r„r[) + p^--^--\r,,r[)] , (8) 

s=±l/2 

whose diagonal part {r[ = ri) represents the one-body density of neutrons or protons. 
We obtain the momentum distributions of protons (t=l/2) or neutrons (i=-l/2) as: 



(9) 



where we have indicated with A/j the number of protons or neutrons. The above definitions imply the following 
normalization of n(k): 

dkn\k) = l. (10) 

We describe doubly closed shell nuclei, with different numbers of proton and neutrons, in a jj coupling scheme. 
The most efficient single particle basis to be used is constructed by a set of single particle wave functions expressed 
as: 

0^y„(rO = Rii,{n)J2 < l^l'^j^ > Yi,,in,)xs{i)xt{i) = Rla,{n)Y^{n{)xt{i) ■ (11) 

In the above expression we have indicated with Y;^ the spherical harmonics, with < | > the Clebsch-Gordan coefficient, 
with Rnij{ri) the radial part of the wave function, and with Y[" the spin spherical harmonics |30 j . 

The uncorrelated OBDMs, those of the IPM, are obtained by substituting in Eq. ([7]) the correlated function 
14" > with the Slater determinant |$ > formed by the single particle wave functions pT|) . We obtain the following 
expressions: 

p*(ri,r;)=^[pr^*(ri,r;)+pr^^*(ri,r;)] , (12) 

S 

pr''\r,,r[) = ^^(2j + l)i?,\y(ri)i?L,(r'i)Pi(cos0n') , (13) 

nlj 
nlj 

In the above equations diy indicates the angle between ri and r[, and Pi and P/ the Lcgendrc polynomials and their 
first derivative respectively. The presence of the second term of Eq. (|12p . the antiparallel spin terms given in Eq. 
is due the jj coupling scheme required to describe heavy nuclei. 
The correlated OBDM is obtained by using the ansatz ([3]) in Eq. ([7]). This calculation is done by using the cluster 
expansion techniques as indicated in 31] and [3^ . where only scalar correlations have been used, and in js^l . where 
the state dependent correlations have been used, but in a Is coupling scheme. We followed the lines of Ref. [33j and, in 
addition, we consider the presence of the antiparallel spin terms and we distinguish proton and neutron contributions. 
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FIG. 2: The difFerence po{ri,r'i) — p{ri,r[), between the proton IPM one-body density matrix of the ^°*Pb nucleus, and that 
obtained with our FHNC/SOC calculations. The two density matrices, have been calculated for Sii/=0. Note that the scale 
here is ten times larger than that of Fig. [1] 



The explicit expression of the OBDM, in terms of FHNC/SOC quantities, such as two-body density distributions, 
vertex corrections, nodal diagrams etc., is given in Appendix [XI The diagonal part of the OBDM is the one-body 
density, normalized to the number of nucleons. Because of this, the momentum distribution satisfies the following 
sum rule: 

Sl = ^Jdkk^ n\k)/T'pHNC = 1 , (15) 

where we have indicated with Tpj^j^^ the kinetic energy of the protons or of the neutrons. We have verified the 
accuracy of our calculations by testing the normalization pop and the exhaustion of the above sum rule for every 
n'(fc) calculated. We found that these quantities are always satisfied at the level of few parts on a thousand in a full 
FHNC/SOC calculation, and even better when only scalar correlations are used. The surface shown in FiglT]represents 
the proton OBDM of the ^ospj-, nucleus, for 6ii'=0. We have shown in [2§| that the SRC lower the one-body proton 
distribution, and also that of the neutrons, in the nuclear center. In order to highlight the effects of the correlations 
on the density matrix, we show in Fig. [5] the quantity Po{'>'i,i"i) — Note that the z-axis scale of Fig. [5] is 

ten times larger than that of FiglTJ It is interesting to notice that the major differences between the OBDMs are not 
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in the diagonal part, but just beside it. The consequences of these, small, differences between the OBDMs on the 
momentum distributions, are shown in Fig. [31 In this figure, we compare the ^^C , ^^O , ^"Ca , ^^Ca and ^°*Pb proton 
momentum distributions calculated in the IPM model, with those obtained by using scalar and operator dependent 
correlations. 




FIG. 3: The proton momentum distributions of the ^^C , ^^O , ''"Ca , ^^^Ca and ^"^^Pb nuclei calculated in the IPM model, by 
using the scalar correlations only (/i) and the full operator dependent correlations (/e). 

The general behavior of the momentum distributions, is very similar for all the nuclei we have considered. Correlated 
and IPM distributions almost coincide in the low momentum region up to a precise value, when they start to deviate. 
The correlated distributions show high momentum tails, which are orders of magnitude larger than the IPM results. 
The value of k where uncorrclatcd and correlated momentum distributions start to deviate, is smaller the heavier is 
the nucleus. It is about 1.9 fm~^ for ^^C , and 1.5 fm~^ for ^"^Pb . We notice that the value of the Fermi momentum 
of symmetric nuclear matter at the saturation density is 1.36 fm~^. 

The results presented in Fig. [3]clearly show that the effects of the scalar correlations are smaller than those obtained 
by including the operator dependent terms. We shall see in the following that this is a common feature of our results. 

Since relatively small differences are compressed in logarithmic scale, we use the linear scale in Fig. |4]to show, as 
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FIG. 4: The proton momentum distributions of the ^'^O and ^"^^^Pb nuclei muhiplied by . The full lines show the IPM results, 
the dotted lines have been obtained by using scalar correlations only, and the dashed lines with the complete correlation. 



examples, the proton momentum distributions for ^^O and ^"^^Pb nuclei, multiplied by fc^. This quantity, multiplied 
by a factor 47r, is the probability of finding a proton with momentum k. We observe that the effects of the SRC on 
the quantity shown in Fig. 0] are basically two. The first one is the already mentioned enhancement at large values 
of k. This effect is less evident here than in Fig. [31 The second effect, hardly visible in Fig. [31 is a reduction of 
the maxima which appear approximately at k—\ fm^^ in both nuclei. These two effects are obviously related, since 
all the momentum distributions are normalized as indicated by Eq. (jlOp . therefore reductions and increases must 
compensate. 

We found that the proton and neutron momentum distributions for nuclei with N ~ Z are very similar. For this 
reason we do not show the neutron momentum distributions of the ^^C , ^^O and '^''Ca nuclei. We compare in the 
panels (a) and (b) of Fig. [Slthe proton and neutron momentum distributions of ''^Ca and ^"^Pb . The thicker lines 
show the results of our FHNC/SOC calculation, the thinner lines the IPM momentum distributions. 

The figure shows that, in our calculations, the differences between protons and neutrons momentum distributions 
are more related to the different single particle structure than to the correlation effects. The main differences in 
the two distributions is in the zone where the n(k) values drops of orders of magnitudes. This corresponds to the 
discontinuity region of the momentum distribution in the infinite systems, which is related to the Fermi momentum. 
In a finite system, the larger number of neutrons implies that the neutron Fermi energy is larger than that of the 
protons, and, consequently, the effective Fermi momentum. For this reason, the neutrons momentum distributions 
drops at larger values of k than the proton distributions. In the panels (c) and (d) of Fig. [5l we show the quantity 

n*(fc)-n*(fc) 

n*(/fc) + n*(fc) ' ^ ' 

where nQ{k) indicates the uncorrelated momentum distribution. This quantity is useful to point out the effects of 
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FIG. 5: In the panels (a) and (b) we show the protons (full lines) and neutrons (dashed lines) momentum distributions of the 
*®Ca and '^"^Pb . The thick lines show the results of our calculations, the thin lines the IPM results. In the panels (c) and (d), 
we show the weighted difference (|16|l between uncorrelated and correlated momentum distributions. As in the upper part of 
the figure, the full lines show the protons results and the dashed lines those of the neutrons. 



correlations. We see that in the low k region A is almost zero. After the discontinuity region A reaches an almost 
constant value around minus one. This behavior indicates that in the low k region the momentum distribution is 
dominated by single particle dynamics. The differences between protons and neutrons at low k are due to different 
single particle wave functions. In the higher k region the correlation plays an important role. We observe that protons 
and neutrons A are very similar, and this indicate that the effect of the SRC is essentially the same for both kinds 
of nucleons. Our results are in agreement with the findings of Ref. 134], where n(fc) of asymmetric nuclear matter is 
presented. There is however a disagreement with the results of Ref. ^35l] . where, always in asymmetric nuclear matter, 
correlations effects between protons were found to be stronger than those between neutrons. 

The increase of the momentum distribution at large k values, induced by the SRC is a well known result in the 
literature, see for example the review of Ref. (36j . The momentum distributions of medium-heavy nuclei, have been 
usually obtained by using approximated descriptions of the cluster expansion, which is instead considered at all orders 
in our treatment. The importance of a complete description of the cluster expansion is exemplified in Fig. [6l where, 
together with our results, we also show the results of Ref. [s^, obtained by truncating the cluster expansion to the 
first order in the correlation lines. In both calculations the same correlation functions and single particle basis, those 
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FIG. 6: Proton momentum distribution of ^''O in various approximations. The thick lines are those of the analogous panel of 
Fig. [3l The thin lines have been obtained by using the first-order expansion method of Ref. [3?! ]. 

of Ref. [29| . have been used. The results obtained with the first order approximation, provide only a qualitative 
description of the correlation effects. They show high-momentum enhancements which, however, underestimate the 
correct results by orders of magnitude. 

III. NATURAL ORBITS 

The natural orbits are defined as those single particle wave functions forming the basis where the OBDM is diagonal: 

nlj 

In the above equation the c^;^ coefficients, called occupation numbers, are real numbers. In the IPM, the natural 
orbits correspond to the mean- field wave functions of Eq. pT|) . and the c^^i^ numbers are 1, for the states below the 
Fermi surface, and for those above it. 

In order to obtain the natural orbits, we found convenient to express the OBDM of Eq. (fT7| as: 

p*(ri,riO = ^*(ri,riOp*(ri,riO + B*(ri,riO , (18) 

where p*(ri,ri/) is the uncorrelated OBDM of Eq. p^ . and the other two quantities are defined as: 

A\v,,v[) = 2C^^soc(ri)Ci,soc(ri)exp[7VL(ri,ri)] + 

2C*(ri)C*(r;)exp[7VL(ri,rl)]^A^-A'=iV*^,p(ri,rl) , (19) 

P>i 

B*(ri,r'i) = -2Cf,,soc(ri)C^f..soc(r'i)exp[7VL(ri,rl)K^^(ri,r'i) - 
2C*(ri)C^(r'i)exp [A^L(ri, r'l)] 

x5]A'=A'={7VL.p(ri,r;X^^(ri,r;)-f<^^^p(ri,r;)} . (20) 
P>i 

The meaning of the ww, Wc^c labels used in the above equations have been defined in (3^ where the index k has been 
defined as p = 2fc + / — 1 with Z = 0, 1 and p — 1, . . . , 6. The detailed expressions of the vertex corrections C and of 
the nodal functions N are given in Appendix VK\ 

We expand the OBDM on a basis of spin spherical harmonics Y[" , Eq. (fTTj) . 

p*(ri,r;) = ^ [A^{rur[)+B\^{r,y,)\ Y*^-^ {n,)Y^ (21) 

Ijm 
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where J7i and ^[ indicate the polar angles identifying ri and r[. The explicit expressions of the A and B coefficients 
are: 



Al^in,r[) - (2^ + 1) ^ (2Z2 + l)(2j2 + l) 



E 



I h h 





h h 3 
I 1/2 h 



with 



and 
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- J dnA'{ri,r[)Piicoseir 



BUri,r[) = J d{cosen')B\n,r[)Pi{cosl 
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) 



(22) 



(23) 



(24) 



In the above equations we have used the 3j and 6j Wigner symbols [30j. The term A depends on both orbital and 
total angular momenta of the single particle, I and j respectively, and the term B depends only on the orbital angular 
momentum I. 

As it has been done in Refs. and [1^ we identify the various natural orbit with a number, a, ordering them 
with respect to the decrea sing value of the occupation probability. The general behavior of our results is analogous 
to that described in Ref. [33| where a system of ^He drops composed by 70 atoms have been studied. The orbits 
corresponding to states below the Fermi level in the IPM picture, have occupation numbers very close to unity for 
a — 1, and very small in all the other cases. As example of our results, we show in Figs. [7] and [8] the proton and 
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FIG. 7: Occupation numbers of the proton natural orbits of the Ca nucleus, having a = 1. The dashed line indicates the 
IPM values. The black bars show the values obtained with the scalar correlation and the gray bars those values obtained with 
the full correlation. 



neutron occupation numbers for the natural orbits witha = 1 of the '*^Ca nucleus. In the figures, the IPM results are 
indicated by the dashed lines. The black bars show the values obtained by using scalar correlations only, the gray 
bars those obtained with the complete operator dependent correlations. 

The correlated occupation numbers are smaller than one for orbits below the Fermi surface, and larger than zero 
for those orbits above the Fermi surface. This effect is enhanced by the operator dependent correlations. We observe 
that, for the states above the Fermi surface, the gray bars are larger than the black ones, indicating that also for these 
states the operator dependent correlations, produce larger effects than the scalar ones. 

We show in Fig. [9] some a — 1 natural orbits for three neutron states in ^^Ca . In this figure, we compare the 
IPM results (full lines) with those obtained with scalar correlation only (dotted lines), and with the full operator 
dependent correlation (dashed lines). The effect of the correlations is a lowering of the peak and a small widening 
of the function. Despite the small effect, it is interesting to point out that this is the only case where we found that 
the inclusion of operator dependent terms diminishes the effect of the scalar correlation. This fact is coherent with 
the results on the density distributions shown in Ref. [2^. In Tab. Uwe show the occupation numbers of the ^^O 
protons and neutrons natural orbits also for a > 1, and we make a direct comparison with the results of Ref. [3^ . 
As already said in the discussion of the ^*Ca results, the inclusion of the state dependent correlations increases the 
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FIG. 8: The same as Fig. [7]for the occupation numbers of the neutron natural orbits of the "^^Ca nucleus. 
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TABLE I: 


Protons (p) and neutrons (n^ 


natural orbits occupation numbers for 


^''O . The PMD values 


are those of Ref. 





differences with respect to the IPM. The occupation numbers of the orbits below the Fermi surface are smaller than 
those obtained with scalar correlations only. The situation is reversed for the orbits with a > 1 or above the Fermi 
level. For the states below the Fermi surface, our full calculations produce correlation effects slightly larger than those 
found in [3§|, whose results are closer to those we obtain with scalar correlations only. For orbits above the IPM 
Fermi surface, our occupation numbers are always smaller than those of Ref. [3^ . 



IV. QUASI-HOLE WAVE FUNCTIONS AND THE SPECTROSCOPIC FACTORS 

The quasi-hole wave function is defined as: 

' [< "fniMA - mmrmiA - 1) X 'fiA)\'f{A) 

where \'^nijm.{^, 1) > ^-nd |^'(1, A) > are the states of the nuclei formed by A— 1 and A nucleons respectively, 

and P4 is the isospin projector. In analogy to the ansatz ([3]), we assume that the state of the nucleus with A ~ 1 
nucleons can be described as: 

^nl.miA - 1) = J^il, A - 1)$„,,™(1, ...,A ~ 1) , (26) 

where ^nijm(XT--T A — 1) is the Slater determinant obtained by removing from $(l,...,v4) a particle characterized 
by the quantum numbers nljm, and the correlation function T is formed, as indicated in Eq. by the two-body 
correlation functions fp obtained by minimizing the A nucleon system. In an uncorrelated system the quasi-hole wave 
functions coincide with the hole mean-field wave functions ()lip . 

We are interested in the radial part of the quasi-hole wave function. We obtain this quantity first by multiplying 
equation (|25p by the vector spherical harmonics Y™(i7), then, by integrating over the angular coordinate fl, and. 



11 



finally, by summing over m. It is useful to rewrite the radial part of the quasi-hole wave function as [s^: 



where we have defined: 



and 
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0.95 


0.87 


0.81 


0.95 


0.83 


0.80 


0.94 


0.83 


0.77 


0.96 


0.90 


0.86 


0.96 


0.90 


0.85 


0.94 


0.84 


0.79 


0.96 


0.92 


0.87 


0.94 


0.92 


0.86 


0.94 


0.86 


0.80 


0.95 


0.90 


0.85 


0.96 


0.90 


0.84 


0.94 


0.84 


0.79 
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0.86 
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0.95 
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0.86 


0.80 
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0.82 
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0.88 


0.83 














0.94 


0.88 


0.82 














0.95 


0.89 


0.83 














0.94 


0.90 


0.86 














0.95 


0.89 


0.83 














0.95 


0.90 


0.85 



TABLE 11: Proton spectroscopic factors of the ^^C , ^""O , '^'^Ca , "^^Ca and ^°^Ph nuclei. We present the results obtained by 
using the scalar correlation only (/i), the first four operator channels of the correlation (/4) and the full correlation operator 

Following the procedure outlined in Ref. [s^, we consider separately the cluster expansions of the two terms Afj^ 
and X^, where we have indicated with a the set of the nljm quantum numbers. We obtain for A"^ the expression: 



^lir) = Cl;"soci^)[RUir)+ d^nRii^{n)Pi {cos 



and for A/"^ the expression: 



+p*'"(r, ri) - iV*;" (r, n) - 7V*'"(r, n) + X'soci^, ri) 



(30) 



d'<'::ai^){mrr+ I dViC(r)0l(ri) 



x2{g**,"(r, ri)C*;;,(ri) [ - p*-"(r, n) + 7V*f (r, r,) 

(r , r 1 ) - (r , r 1 ) - N^^; (r, ri ) + AA* (r , r i ) } ) , 



(31) 
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12^ 






16q 




/l 


U 


h 


fi 


.A 


/6 


0.96 


0.95 


0.91 


0.95 


0.90 


0.85 


0.96 


0.96 


0.94 


0.96 


0.93 


0.89 








0.96 


0.93 


0.89 





«Ca 
















h 


/4 


/6 


/i 


/4 


/6 


/i 


/4 


/6 


0.93 


0.84 


0.78 


0.93 


0.86 


0.80 


0.92 


0.85 


0.80 


0.95 


0.87 


0.82 


0.94 


0.88 


0.83 


0.93 


0.85 


0.80 


0.95 


0.87 


0.81 


0.94 


0.88 


0.82 


0.93 


0.85 


0.80 


0.96 


0.90 


0.86 


0.95 


0.90 


0.86 


0.93 


0.86 


0.82 


0.96 


0.92 


0.87 


0.95 


0.92 


0.87 


0.93 


0.88 


0.84 


0.95 


0.90 


0.85 


0.95 


0.90 


0.86 


0.93 


0.86 


0.82 
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0.93 
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0.95 


0.94 
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0.94 


0.93 


0.90 














0.95 


0.94 


0.90 



TABLE III; The same as Tab. Ullfor neutron states. 



The expressions of the functions AfgQfj{r,ri), XgQ(j{r,ri), are 



Xpoc(.r,ri) ^Y."^' E [(l-'5fc,i)'^2V-i,2fe-i(r,ri) 

fc=l t2=p,n 
k=l t2=p,n 

where the indexes t refer to the isospin, and wc have defined: 

AA^Vnr,ri) = [/»*^;-"(r,r05^i,*^'"(r,rOC*-"(rO(-p*-"(r,ri) + A^*-"(r,ri)) 



^k2 



(32) 



(33) 



(34) 



(35) 



where k2 — 1,2,3 for q = 1,3,5. The expressions of the other terms are given in Appendix [XI All the quantities used 
in the above expressions depend on the set of quantum numbers a characterizing the quasi-hole state, since we have 
written the various equations by using [40| : 



p*'"(r,ri) = p*(r,ri)-(/.l*(r)(/.l(ri) . 
The knowledge of the quasi-hole functions allows us to calculate the spectroscopic factors: 

5*y = / cirr2|V'^y(r)|2 . 



(36) 



(37) 
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FIG. 9: Natural orbits for some neutron states in Ca . The full lines show the IPM orbits, the dotted lines those obtained 
with scalar correlations only and the dashed lines those obtained with the complete operator dependent correlation. 

The proton and neutron spectroscopic factors for the ^^C , ^^O , '*°Ca , '^^Ca and ^osp]-, jmdgj are given in Tabs. 
im and mil for each single hole state. In these tables we compare the results obtained by using scalar correlations 
(/i), with those obtained with the four central channels (fi) and with the full correlation (/g). The inclusion of the 
correlations produce spectroscopic factors smaller than one, the mean-field value. The /g results are smaller than 
those of /4, which are smaller than those obtained with /i. 

We found that the effect of the correlations becomes larger the more bound is the state. This fact emerges by 
observing that for a fixed set of Ij quantum numbers the spectroscopic factors increase with n, and, at the same time, 
that the values of the spectroscopic factors become larger when n and the Ij values increase. 

The values of the spectroscopic factors depend on the choice of the single particle basis. As we have already said 
in the introduction, our results have been obtained by using the Woods-Saxon single particle bases given in Ref. [1^ . 
This basis has been chosen in order to reproduce the single particle energies around the Fermi surface and the charge 
distribution of each nucleus considered. The correlation function has been fixed by the minimization procedure 
We tested the sensitivity of our results to different single particle basis, by calculating ^^O and ^°Ca spectroscopic 
factors by using with the Harmonic Oscillator and Woods-Saxon single particle wave functions of Ref. [41| , fixed by 
a global minimization of the energy. Despite the remarkable differences between the various single particle basis, we 
found that the maximum variations in the values of the spectroscopic factors is of about the 5%. This value is smaller 
than the variations produced by the different terms of the correlations, shown in Tabs. HIl and Hill This indicates that 
our results are more sensitive to the SRC than to the choice of the single particle basis. 

As example of correlation effects on the quasi-hole wave functions, we show in Fig. [TO] the squares of the proton 
Ssi/2 and neutron ipi/2 quasi-hole wave functions for the ^"^^Pb nucleus. The correlations lower the wave function in 
the nuclear interior. Also in this case, the effect of the correlations increases together with the number of operator 
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FIG. 10; Proton 3si/2 and neutron 3pi/2 quasi-hole functions, squared, of the '^"^Pb nucleus. The various lines show the results 
obtained by using different type of correlations. 

channels considered. 

In Fig. [Tl]we show with a gray band the difference between the empirical charge distributions of ^°^Pb and ^osxi 
[l9l |. The dashed dotted line, labeled as IPM, has been obtained by considering that the difference between the two 
charge distributions can be described as a single 3si/2 proton hole in the core of the lead nucleus. This curve has been 
obtained by folding the IPM result of Fig. [TO]with the electric proton form factor in its dipole form. In a slightly more 
elaborated picture, the ground state of the '^'^^T\ is composed by the 3si/2 proton hole in the •^"^Pb ground state, plus 
the coupling of the 2d^/2 and 2^3/2 proton levels with the first 2+ excited state of 206p|-, [42^|4^. This description of 
the 205^1 charge distribution, shown by the dotted line in the figure, is still within the IPM framework. The dashed 
line has been obtained by adding to the dotted line the core polarization effects produced by long-range correlations. 
These effects have been calculated by following the Random Phase Approximation approach of Refs. [3,|4^. The 
full line has been obtained when our SRC effects are also included. 

The various effects presented in Fig. Illl have been obtained in different theoretical frameworks, and the final result 
does not have any pretense of being a well grounded and coherent description of the empirical charge differences. The 
point we want to make by showing this figure is that the effects of the SRC are of the same order of magnitude of 
those commonly considered in traditional nuclear structure calculations. 

V. SUMMARY AND CONCLUSIONS 

In this work we have extended the FHNC/SOC scheme in order to calculate the OBDM's, the natural orbits and 
the quasi-hole wave functions of finite nuclear systems non saturated in isospin, and in jj coupling representation 
of the single particle wave function basis. Our results have been obtained by using the many-body wave functions 
obtained by minimizing the nuclear hamiltonian with the two-body realistic interaction Argonne Ug and the associated 
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FIG. 11; Differences between charge density distributions of ^"®Pb and ^'^^Tl. See the text for the explanation of the various 
Unes. 

thrcc-body interaction Urbana IX. The calculations have been done by using operator dependent correlations which 
include terms up to the tensor ones. 

We found that the correlations enhance by orders of magnitude the high-energy tail of the nucleon momentum 
distribution. The occupation numbers of the natural orbits below the Fermi level, are depleted, and the opposite 
happens for those above the Fermi level. Also the values of the spectroscopic factors are depleted with respect to the 
1PM. A reliable comparison between our spectroscopic factors with the empirical ones requires the description of the 
reactions used to extract them, and this is part of our future projects. 

We have shown that the results of models considering expansions up to the first order correlation lines, provide 
only qualitative descriptions of the SRC effects. In the description of the charge density difference between 206pj^ 
and ^'^^Tl, the SRC effects are of comparable size of those commonly considered in traditional nuclear structure 
calculations based on effective theories. 

A general outcome of our study, is that the effects of the correlations increase with the complexity of the correlation 
function. This means that operator dependent correlations enhance the effects produced by the scalar correlations. 
This not obvious result, is valid in general, not always. We have shown in Ref. that scalar and operator 

dependent correlations have destructive interference effects on the density distributions. We found in the present 
study an analogous behavior regarding the natural orbits. These quantities are related to the density distributions. 
It seems that the effects of the SRC are rather straightforward on quantities which involve two-nucleons, while they 
are more difficult to predict on quantities related to single nucleon dynamics. On these last quantities, however, these 
SRC effects are very small, usually negligible. 

In our calculations the nuclear interaction acts only in defining the many-body wave functions by means of the 
variational principle ^ , more specifically, in selecting the correlation function It is therefore difficult to disentangle 
the role played by the various parts of the interaction, e.g. the tensor part of the three-body force, on the quantities 
we have studied in this article. We have instead evaluated the effects of the various parts of the correlation function. 

In this work, we have highlighted a set of effects that cannot be described by mean field based effective theories. 
The description of the nucleus in kinematics regimes where these effects are relevant, requires the use of microscopic 
theories. 
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APPENDIX A 



For sake of completeness, we give in this appendix the detailed expression of the OBDM for finite nuclear systems 
not saturated in isospin, and in jj coupling scheme of the single particle wave function basis (jlip . The notation for 
the nodal functions N and for the vertex corrections C is that used in Ref. ^SSj. The indexes ti, t2, indicate protons 
and neutrons, and the subscript j is related to the antiparallel spin states. 

For the correlated OBDM we obtain the expression: 



p*(ri,riO = 2C*,soc(ri)C*,5oc(ri')e~--(-^'"-^'^ 



p*(ri,riO -iV* (ri,riO 



(Al) 



P>i 



In the above equation, k has been defined as in Eq. ((20l) . and we have used A'^"^^^'^ = 1 — ^^^3, and A'^~^''^'^ = 1, 3, 6. 
In the following we shall calculate the expectation value of the isospin operator sequence: 



by considering that 



with 



a„+i = 3(1 - a„) 
By using the above equations we have that: 



{Ti ■ Tj)" = a„ + (1 - a„)Ti • Tj , 

ao = 1 . 



and 



Xo'*' = l , x^*' = 2<5t,t, ~ 1 and = 2a„ - 1 + 2(1 - a.^'^ut. 



The expressions of the vertex corrections are: 

C^oc(ri) = C*Hri)fl + C/::soc(ri) 



fe=l t2=p,n 



where we have defined 



h{r 12) 

The expressions of the two-body distribution functions for p > 1 are: 



ffL,p(i"i,ri') 
9ttp(^i,r2) 

9l%(ji.^2) 
5!;ec,p(i'i,r2) 



3^(1-1, i-i')^c^,L,p(i-i, 1-1') , 
/i!;,pHri,r2)£f (ri,r2) 

<^i,p(ri,r2)+XX(ri,r2), 

(r 1 , r2)g'f^ (ri , r2 ) + (r 1 , V2)Nl% (n , r2 ) 

^L\p^ (ri, r2)£^^,(ri, r2) + (n, r2)7V*^^,,p(ri, ra) 



(A2) 



(A3) 



(A4) 
(A5) 

(A6) 
(A7) 
(A8) 
(A9) 
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Finally the nodals functions are expressed as: 

2) - <^„.(,),,(1, 2) + 2) + 2) + Nl]^{l, 2) 



(AlO) 



with m,n = c,Wc. The separation of the above nodal diagrams in four terms, corresponds to the classification in the 
XX, xp, px and pp parts [29|, and it has been applied to the quantities -^^^j) pgri^' 2) defined in the following. 

3 

^mm, 2fci-l(-'-' 1 ) = 5Z 5Z [^mm, 2fci-l, 2fc2-l, 2fc3-l(-'-' 1 ) 
^2^3 = 1 *3=P," 

+ Xl^ ^ [^mm,2fci-l,2fc2,2fc3-l("'"' ^ ^mm, 2fci -l,2fc2 - 1, 2/03 ' ^ )] 
3 

-^rnin,2fci (1' 1 ) = ^ ^ ^2^ ^rnm,2fci ,2fc2 ,2fc3 ' ^ ) ' 



fc2,fc3 = l *3=P," 
3 

-^wrif2fci-l(-'-' 2) = ^ ^ [^wn,2fci-l,2fc2-l,2fe3-l(-'-' 2) 
fc2fc3 = l *3=P," 

I,2fc2,2fc3-1 I,2fc2-l,2fe3 

(1,2) 
3 

^wn,2fei(-'-' 2) = ^ ^ -^wn,2fci,2/c2,2fc3(-'-' 2) , 
^2,^3 = 1 *3=P," 
3 

^c(j),2fci-l(-'-' 2) = ^ ^ [^Jec(i),2fci-l,2fe2-l,2fe3-l'-"'"' 2) 

fc2,fc3 = l *3=P," 

+ Xl^ ^ [^JecO"),2fci-l,2A:2,2fc3-l(-'-' 2) + ^J^c0),2fei-l,2fe2-l,2fc3 (1' 2)] 
3 

^(!.ecO%2fci(l'2) X! X! ^i^'cc'(j),2fei,2fe2,2fc3'^"'"' 2) • 

fc2,fc3 = l *3=P," 



with m = u;, Wc and n — d,e. 



A^*tp.^.(ri,r2) = [^*t,(^i''-3)^il2^'^^^!,.(r3)|X*?*Ur3,r2) 

'■3)^il2^''^e?,.(r3)|^*^:Mr3, r2) + A^etMr3, rs)" 



also in the above equations we used n = d,e. In the following equations we have that m,n = c, w, 

r A*^^ 

(1 - 5r^i) X 

( ^ ^ ^ /\fe2 J ^ ^ ^ ■ 

-'^mc(l, 3)^l|2^ —^Cl'^qr{3)\Xlf^ j,{3, 2) + A^'^^^^rl^, 2) + iV*^ ,,(3, 2) 



4,i42i43i{ 3)ic*V(3)|^*^,(3, 2) + Nll^^{3, 2) + TV^^/S, 2)^ 

+ (1 - [X^,,(l, 3)ic*V(3)|X*^^. ,(3, 2) + 7VX-,r(3, 2) + A^^i^,, .(3, 2)] } , 



(All) 
(A12) 

(A13) 

(A14) 
(A15) 

(A16) 
(A17) 



(A18) 



(A19) 



(A20) 
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r A'^^ 
2) = X'^,Jl,3)e^l^^''^ -^Cl%,m - (3, 2) + Nlia, 2) + N'; (3, 2) 



<p?,.(l,2) 



Ar*i*3 . (19) 



r A'^^ 
+ (1 - Sr.i) [X^,(l, 3)^1^3^'=^ — C*V(3)l^^c?,.(3, 2) + 7V;^,(3, 2) 



^fc,i'5;c.i^fc3i{ 3)-C:;',,(3)| - p*}(3, 2) + iV,%.(3, 2) + iV;](3, 2) 

+ (1 - ^.,1 ) [xl/1 , 3) ic*;^,, (3) |A^*^,- , (3,2) + Nil, (3, 2)] } , 

Po^ (1, 3)C^3t'^ —ClU^cU^' 2)_ 
Po^(l,3)fi1^'-^— (C*V(3) ~ l)|iV*;^.(3,2) - 6r,ipl'{3,2) 
+Sk,iSk,iSk,i{ [pI){1, 3)^C*f,,(3)|iV*^^. ,(3, 2) 
+ [pojih5)l{Cl:,A3) - 1)\n;1^{3,2) - 5.,ip*}(3,2)] } , 
4a'5fe2i43i{ 3) (3) 1^*^,(3, 2) + iV*f,,^.(3, 2) + iV;?,^.(3, 2) 

(1 - Sr,i) [X'„UI, 3)ic*f,,(3)|X*;^^. ,(3, 2) + 2) + iY^-^,. ,(3, 2)" 

'X'^cj.,ih 3)\cl^grmKU3, 2) + NlU3, 2) + 7V*?,(3, 2)] + 
(1 - Sr,i) [xL/l, i)lcl^grmxXr{3, 2) + A^*;^,,.(3, 2) + iV;,^,,,(3, 2)] } 



NXlp,rih 2) = 4a42i43i{ [^^,,(1, 3)ic*f,,,(3)| - p%{3, 2) + iV*^^.(3, 2) + Nl]{i, 2) 
+ (1 - J,,i)k*;e(l,3)ic*;>^,(3)|iV*^^.,(3,2) +iV;^V(3.2) 



+ 



mcj,q 



(1, 3)-C:f,,(3)| - p*H3, 2) + Nl;{3, 2) + iV;n3, 2) 



(1 - Sr.i) [^^c,(l, 3)-C*3^,(3)|7V*^_,(3, 2) + iV*;^,(3, 2)] } 
-4a42i43i { [Po^ (1 , 3) ic^f,, (3) |A^c%-,,. (3, 2)" 
+ [p*^ (1, 3)i(C*f,,(3) - l)|Af;,V(3. 2) - 5r,iP%{i, 2) 
Poi(l,3)ic*V(3)|7V*^,,(3,2) 

PoMl'3)^(C*;^,.(3) - 1)|7V;3,,(3,2) -5.4p*^(3,2)]} . 



(A21) 
(A22) 



(A23) 



(A24) 



(A25) 



The values of the S.ijil^'''^^ coefheients are given in Ref. j46 
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